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ABSTRACT 

The magnetohydrodynamic dynamo equation is derived within general relativity, us- 
ing the covariant 1 + 3 approach, for a plasma with finite electric conductivity. This 
formalism allows for a clear division and interpretation of plasma and gravitational 
effects, and we have not restricted to a particular spacetime geometry. The results 
should be of interest in astrophysics and cosmology, and the formulation is well suited 
to gauge invariant perturbation theory. Moreover, the dynamo equation is presented in 
some specific limits. In particular, we consider the interaction of gravitational waves 
with magnetic fields, and present results for the evolution of the linearly growing 
electromagnetic induction field, as well as the diffusive damping of these fields. 
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1 INTRODUCTION 

Magnetic fields play a vital role in many astrophysical systems and i n cosmology, and the possible sources for the fields on 
different scales has been the focus of immense research o ver the years jKronbere)ll994l : ban fc Wielebinskilbooi |Parkerlll979l : 
IZeldovich. Ruzmaikin fc SokolofJll983l:lBeck et al.lll996|) . In particular, the origin of large scale cosmol ogical fields is still not 
clear l y understood. Although there are well understood mechanisms for amplifying a given seed field JGrasso fc Rubinstein! 
1200ft TwidrowllioO^ . the explanation of the seed sources is more cumbersome, especially in the case of an inflationary sce- 
nario. A wide range of possible physical explanations for estab l ishing a seed field o f high enough field strength has been pre 
posed in the literature (see, e.g . . iHarrisonl lll97ff):lHogan Jl983h: |Vachasp atil Jl99ll):IOaspermi. Giovanini fc Ve nezianol lll99fil 



posed m tiic literature (sec, e.g., Harrison (19/0); iiogan (1983); Vacnasoati (1991); uaspcrim. uiovamm & Vencziano (199o); 
ISid. Jedamzik fc Olintol lll997l):ljovce fc ShaposhnikovHl997ypavis et al.N200ll) : ICalzetta fc Kandud j2002MDimopoulos et ail 



j2002l) : lBetschart. Dunsbv fc Marklundl J2004h : iMatarrese et alJ (120041)). The main mechanism be hind the amplification of the 
seed fields is the dynamo JParkerlll979t IZeldovich. Ruzmaikin fc Sokolofllll983l : iBeck et al.lll99t3) . for which the fluid motion 
manages to amplify even very weak fields on very short time scales. The dynamo mechanism derives from the magnetohydrody- 
na mic (MHD ) appro x imation, and it has been ana l ysed v i gorously within, e.g., expandi n g cosmological models (for a selection, 
see iFermelivl lll980l): Isil. Baneriee fc Chatteriee! Jl99d): iGailis. Frankel fc DettmanrJ jl995l) : Ijedamzik. Katafinic fc Olintol 
il998l) : lBnmdenbmg^Emm references therein). 

Cosmic magnetic fields are observed on all but the largest scales. On the other hand, gravitational waves have yet to 
be directly detected. A lthough strong evidence exists in their favour (such as the PSR 1913+16 binary pulsar observations 
ijHulse fc Tavlorlll97iTh ). the very weak interaction of gravitons with light and matter makes an earth based detection diffi- 
cult. Currently, gravitational wave detectors, both of interferometer and bar type, are coming on-line, and are expected to 
start collecting important astronomical data, opening up a new window for observational astrophysics and cosmology. As 
gravitational wave information is gathered, the detailed study of typical wave signatures will be of major importance. Thus, 
knowle dge of both direct gra vitational wave signature, and indirect signatures, e.g., electromagnetic waves induced by gravity 
waves ((Clarkson et alJboO^h . will facilitate the study of gravitational waves, and is likely to bring new understanding to the 
observations. 

The interaction of gravitational waves, and other general relativistic gravity effects, with electromagnetic fields and plas- 
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(200( j):lTsagasl (2001)). Here, we will derive the general form of the dynamo equation, using the 1 + 3 covariant formalism 
l Enrs&^anElstin998l) . and analyse its coupling to matter and spacetime geometry. This formalism allows for a clear cut 
interpretation of gravitational vs. purely MHD effects, and also presents the influence of the gravitational kinematics on the 
MHD modes for straightforward comparison. Specifically, we will consider the propagation of gravitational waves through a 
resistive magnetohydrodynamic plasma, and derive covariant equations for the propagation of a resulting magnetohydrody- 
namic wave. Here, a word of caution concerning the use of the phrase 'dynamo' is in place. While speaking of the dynamo 
equation, a general evolution equation for the magnetic field in diffusive media is assumed, while the dynamo mechanism is a 
particular application of this general equation. Indeed, since the dynamo mechanism is without question the most prominent 
feature of the general dynamo equation, we would like to stress that although we derive the general dynamo equation, we 
will not apply this to the dynamo mechanism. Instead, we will investigate what could appropriately be termed induction 
fields. The waves presented here represents linearly growing induction fields, and are thus not due to the proper dynamo 
mechanism which instead gives rise to rapid amplification of seed magnetic field due to the vortical fluid motion. Furthermore, 
the effects of a finite conductivity is presented. The equations are analysed numerically, and we discuss possible applications 
to astrophysical and cosmological systems. 



2 BASIC EQUATIONS 

When analysing low frequency magnetised plasma phenomena, MHD gives an accurate and computationally economical 
description. Specifically, a simple plasma model is obtained if the characteristic MHD time scale is much longer than both 
the plasma oscillation and plasma particle collision time scales, and the characteristic MHD length scale is much longer 
than the plasma Debye length and the gyro radius. These assumptions will make it possible to describe a two-component 
plasm a in terms of a one-flu i d description, as the electrons can be c onsidered inertialess (see also, e.g. , |s_ grvin fc Brodinl 
ll2003l) : lMoortgat fe Kuiipersl feooal |2004 : IPapadopoulos et al l (l200ll) : |Papadopoulos. Vlahos fe EspositoJ hooi) ). The one- 
fluid description means a tremendous computational simplification, especially for complicated geometries. Moreover, if the 
mean fluid velocity, the mean particle velocity, and the Alfven velocity is much smaller than the velocity of light in vacuum, 
the description becomes non-relativistic and simplifies further. 



2.1 Covariant theory 

The 1+3 covariant approach relies on the introduction of a timelike vector field u a : u a u a = —1, with which all tensorial objects 
are split into their invarian t 'time' and 'space' parts: scalars, 3-vectors, and projected, symmetric, trace-free (PSTF) tensors 
(see lEllis fc" van Elstl 11199a) for a comprehensive review) . This allows us to write the coupled Einstein-Maxwell equations in a 
relatively simple and intuitive fashion. Here, we introduce consistent 3-vector notation of Euclidean vector calculus, making 
the fully general relativistic equations relatively easy to read. 

In an arbitrary curved spacetime, in units such that c = 1 and G = I/Stt, Maxwell's equations take the form llTsagas fc B arrow! 
ll997HEllis fc van Elstlll993 iMarklund et ahlEooll 

B {a) + curl£ a = - (ZQh ab - a ab + e a6t W) B b - e aba u b E c , (1) 

-£ <Q> +curlB a = Mo j <a> + (%Qh ab - a ab + e abc u c ) E b - e abc u b B c , (2) 
B a E a = 2u> a B a , (3) 
D a B a = -2u> a E a , (4) 

where curl X a = e abc D b X c and X' a ' = h ab X b for any spatial vector X a , and h ab = g a t + u a u b is the metric on the local rest 
space orthogonal to the observer four- velocity u a . Moreover, we have assumed quasi-neutrality, i.e., the total charge density 
satisfies p ~ since the electrons and ions follow the same motion. Furthermore, we will assume that Ohm's law in the form 

ji"-) = _L ( E a +e abc v b B c ) , (5) 



Mo A 

holds, where = h ab j b is the three-current, (/ioA) -1 is the electric conductivity, A is the magnetic diffusivity, and v a is the 
plasma three-velocity. We note that high conductivity implies small magnetic di ffusivity. The case of non-id eal MH D in curved 
space times has been treated sparsely in the literature (with some exceptions jFennellvlll980t Ijedamzik. Katalinic fc Olintol 
1998)), and the possible effects of a finite conductivity due to gravity-electromagnetic field or fluid couplings has therefore 
largely gone unnoticed. 

Thus, taking the curl of the curl of B a and using the commutator relation 

curl (curlB) a = —~D 2 B a + D a (D 6 B b ) + 2e abc B (b) iu c + R ab B b , (6) 
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where R a b is defined by the relation 

R bd = | {£ + A - |6 2 - a 2 + to 2 ) h bd + Il bd - <r <M ) + D (() tt d ) 

-|6(cr6d -Ubd) - u {b u d ) - 2 ^(J (i) fc cr d ) fc +U( b u) d ) \ - 2a k[b uj k d] , (7) 

where S = Skmd + \{E 2 + B 2 ) is the total energy density, and IT 6 = IIfl* id - (E {a E b) + B {a B b) ) is the total anisotropic 
pressure. Taking the curl of Eq. using Eq. Q and @ to replace E a by B a and j^, we obtain 

curl (curlB)" = e abc e cde T> b {v d B e ) ~ (§9 + A -1 ) \B {a) + ^Qh ab - a ab + e abc Lo c ) B b - Xfi e ahc u b j (c) + e abc e cde u b v d B e ] 

-|^oAe abc j( 6) D c e + ^e abc e bde v d B e D c Q - e abc Dj, \a cd ^/i Aj (d> - e def v e B^ - e abc D 6 [e cde ^ d (/ioAj <e) - e efB VfB g ^ 

-e abc D b (e cdc u d B e ) + curl E {a} . (8) 



Furthermore, the last term in Eq. 10, due to the displacement current E^ a \ may be written 

curl£ (a> = -B {a) +E a (9) 

where 

2 r\uab • (ab) , abc \ 7-1 / r\uab ab , abc \ 7-1 1 /~* / <t& 1 abc \ 7-1 

Bra — cr '+€ w/ c \ ijj - on (T +e u c JJm -^B |9fi — cr +e w c ) B b 



—e abc u b |^curli? c - fioj(c) + e cd cU d B e - (^Qh cd — o cd + e cde uj e ) S d j + \e abc E b q c + [2u^ tt a>^ + D^uj^ - (curla) ab j E b 

-e abc u {b ) E c + e abc (a bd + e bdB UJ e ) D d £ c (10) 

by commuting spatial and time like covariant derivatives and using Maxwell's equations. The commutation of derivatives 
introduces curvature effects into the expression. Here q a — q^ uiA + e abc EbB c is the energy flux due to fluid and electromagnetic 
(i.e., Poynting flux) contributions. The latter derives from the gravitational self-interaction of the electromagnetic field, and 
can in many cases safely be neglected. We note that all terms in H a defined by Eq. (HUH are curvature contributions. The 
displacement current is normally neglected in the MHD approximation, since it corresponds to high frequency phenomena. 
Here we see that the presence of gravity alters this interpretation of the displacement current. We will be interested in low 
frequency phenomena, and we therefore assume curli5' a ' ~ H a . 

Thus, equating @ with the expression @, we obtain the general relativistic dynamo equation (see also iTsagasI l|2004h 
for the source free general relativistic electromagnetic wave equations) 

g(a) _ e abc ecdeDb (i, d B e ) - \D 2 B a = -§AeB <a> + 2Xe abc io b B (c) + 2AD a [w t ^ Xj {b) - e bcd v c B d ^ - XR ba B b 

f -t i 2 \ s\\ \ ( 2 ai ctb ab , abc \ n \ abc • ■ , abc ■ d 7->e 

-(1+ jAB) I jB/i — a +e lo c j Bb — fJ-oXe u b j( c ) + e e cd£ u b v B \ 
~|/ioA 2 e ai,c j( b )D c 6 + ^Xe abc e bd cV d B e DcQ - \e abc ~D b \a cd (jioXj^ — e def v e Bf^j j 

-Ae a6c D b \e cde u A [^\j {e) - e ef3 v f B g )] - Xe abc D b (e cde u d B e ) + AH a (11) 

where we have used Maxwell's equation Q and Ohm's law ©. The terms on the right hand of the dynamo equation il l 1 1 
represents the influence of general relativistic gravity, while the left hand side gives the normal dynamo action from the fluid 
vorticity e abc D b v c . 

Sometimes it may be advantageous to keep the electric field in the gravitational right hand side of Eq. I I U , in which case 
we have 

£<a> _ e abc eede D b (v d B e ) - \F> 2 B a = -|AeB <a> + 2Xe abc io b B (c) + 2AD a (u b E b ) - XR ba B b 
- (1 + |A6) [ ( K l®h ab - o ab + e abc Lo c ) B b - e^iibEc] - le abc E b D c Q - Xe abc D b (a cd E d ^ 

-Xe abc D b (e cde uj d E e ) - Xe abc D b (e cd cU d B e ) + XE a . (12) 



2.2 Three-dimensional vector notation 

The general relativistic dynamo equation may also be formulated in terms of three-dimensional vector notation. For every 
spacelike vector X a or PSTF tensor Y ab we denote the corresponding three-dimensional vector or tensor according to X and 
Y. With these definitions, Eq. 112H takes the form 

B-VxpxB) - AV 2 B = -|A6B + 2A Co x B + 2AV(£ • E) — XB • R — (l + §A6) (%®B -B-B + BxCu-axE^J 
-|A£ x Ve - AV x (a ■ E) - AV x (Co x E) - AV x (a x B) + AH, (13) 
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where we have introduced the three-dimensional acceleration vector a corresponding to ii a , and V corresponds to D a . Dot 
and cross products are defined in the obvious way. Furthermore, the terms due to the displacement current takes the form 

a = - (§ QB -a ■ B + B X uj - (fO-B -a ■ B + B x u5\ - 'a x E + §9V x E + ■ V) x E - (Co x V) x E 

-ax [VXB - ^ J+ ax B — (^@E - a ■ E + E x uA] + \E x q- H ■ E, (14) 

where we have introduced the magnetic part of the Weyl tensor H ab = —2u^ a io b ' > — D l * a Lo b ' > + (curia)" 6 , which signifies the 
presence of gravitational waves or frame dragging effects. 

2.3 The equations of motion 

The dynamo equation contains the centre of mass fluid velocity v a = {S^ e )V? e \ + <%)«(;) )/(<£(e) + <%)) an d the current 
j' a ' = p( e )f" e ) + p^vf^, e (i) denoting the electrons (ions). Using quasi-neutrality, the evolution of the total fluid energy 
density (dropping the index denoting the fluid) S = Su,\ + $U\ and the centre of mass velocity is given by 

S + D a (£v a ) = -(6 + ii a v a )S, (15) 

and 

S (i {a) + v b D b v a ) = - (l9« a +u a + o\v b + e abc u; b v^ S + e abc j b B c , (16) 

respectively, in the cold plasma limit. Moreover, the current j a can be expressed in terms of the magnetic field via Maxwell's 
equation and Ohm's law JjJ. 

We note that the Eq s. 115[ and l|16fl can easily be generalised to incorporate anisotropic and/or viscous effects (see, e.g., 
ISubramanian fc Barrowl (|l99St) '). 

3 SPECIAL CASES 

In order to extract some information from the rather complicated relativistic dynamo equation, we shall present some special 
cases. 

3.1 Infinite conductivity 

The large bulk of literature on general relativistic MHD has treated the case of ideal MHD, i.e., infinite conductivity. As 
shown below, the general dynamo equation greatly simplifies when this assumption is made. 

As a — > oo, the diffusion of the magnetic field lines decreases, i.e., A — > 0, and we are left with 

ri(a) abc j-\ t dj->e\ / 2nLnf' o,b . abc \ 7-, abc • dne fir>\ 

B w — e e c deUb(v B ) = — I gBft —a + e uj c J B b - e e c deU b v B , (17) 

or 

B-Vxpx5) = - -B ■ B + B x Cu + ax (tf x B\ 1 , (18) 

which of course also can be found directly from Maxwell's equations and the vanishing of the Lorentz force. Infinite conductivity 
is a common assumption in cosmology, for example. 

3.2 Effects of expansion/collapse 

As in the case of infinite conductivity, the case of curved spacetime MHD has largely resorted to analysing the effects of 
expansion/collapse, at which the gravitational-electromagnetic field coupling is kept at a minimum. It also gives a good 
description of some fundamental aspects of cosmological MHD phenomena. 

The simplest case where the nontrivial effects of expansion occur is in homogeneous and isotropic fluid background. The 
homogeneity and isotropy significantly simplifies the Einstein equations, and the spacetime is described in terms of the scalar 
quantities S ', p, O, and A. Therefore, equations II U and 1131 simplify to 

B a - e abc e cde D b {v d B e ) - XD 2 B a = -§AeB a - XR ba B b - (l + §A6) f QB a + \E a , (19) 

and 

B - V x (v x B) - \V 2 B = -|A6B — XB • R — (l + §A9) |6B + AS, (20) 
respectively, where R a ;, = §(<? + A — ^Q 2 )h ab is the Ricci three-curvature. The effects of the displacement current becomes 

E a = -§OB a - QB a - l<d 2 B a , (21) 
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i.e., 

H= -l&B-BB- |e 2 B (22) 

Here we have neglected the gravitational self-interaction of the electromagnetic field. Thus, from Eq. 1191 and 121)1 . together 
with 112 U and 12211 . we obtain 

B a - e ahc e cde D b (v d B e ) - XD 2 B a = ~\Q{2B a + 5\B a ) - § A [S - 3p + 4A + §9 2 ) B a , (23) 

and 

B — V x (v x B) — AV 2 B = -|6(2B + 5AB) - |A {£ - 3p + 4A + §9 2 ) B, (24) 
respectively. Thus, we see that for collapsing solutions (O < 0), there will be a growing magnetic field. 



3.3 The effects of rotation 

In a rotating space time, i.e., u a 7^ 0, Einstein's equations shows that the vorticity has as its source the acceleration u a . 
Godel's universe does not require this, as it is spacetime homogeneous, but for the sake of completeness we will keep both 
spacetime acceleration and vorticity, while neglecting the expansion and shear (thus only considering only rigid rotation), in 
the equations presented below. They read 

B {a} - e abc e cde D b (v d B e ) - XD 2 B a = 2Xe abc cj b B {c) + 2AD a [u b (/W b> - e bcd v c B d ^\ - \R ba B b - Xe abc D b (e cde u d B e ^ 

I abc t~> \ abc • ■ , abc ■ d 7->e \ abc-w-\ d I \-(e) efq n \ 1 \ 1 — <& /or\ 

— — e oj b B c — /ioAe u b j^ + e e cde u b v B — Ae U b \e cde uj \^ Q Xy ' - e ja v f B g ) \ + \=: (25) 

or 

B-Vx(txB) - AV 2 B = 2Aux B + 2AV [dJ • [fioXj-vX s)J - XB ■ R - \S x u - a x (fi \j - v X Sj ] 

— AV x \io x (^/ioAj - « x fl)] - AV x (a x B^j + AH. (26) 

We note that the contribution H a from the displacement current is only slightly simplified compared to the generic case. 
Thus, as in many problems of general relativity, rotation of spacetime gives rise to highly complex equations, and the dynamo 
equation is no exception. 



4 GRAVITATIONAL WAVES 

Linear gravitational waves on a homogeneous background can be covariantly defined as tensor perturbations satisfying D a (r ab = 
~D a E ab = D a H ab = 0, thus being transverse and traceless. In order to clearly see the effects of gravitational waves on the 
induction of electromagnetic fields, we consider this perturbation on a Minkowski background. In this case, the dynamo 
equation becomes 

B (») _ e abc e cde D b (v d B e ) ~ XD 2 B a = -\R ba B b + a ab B b - \e abc D b (a cd E d ) + \E a , (27) 

or 

B - V X (v X B) - AV 2 B = -XB ■ R + a ■ B - AV x {a ■ E) + AH, (28) 

where we have neglected all products of gravitational variables. 
The contribution from H a becomes 

E a = &(ab ) Bb + ^b^^ + ^bc^jyd^ (curl£T) a^ 6 (2g) 

Equation J7J gives an expression for the Ricci three-curvature in terms of the shear of the perturbation. Thus, using 12911 . 
Eqs. and give us 

B {a) - e abc e cde D b (v d B e ) - XD 2 B a = a ab B b + X \2d {ab) B b - e abc D b (a cd E d ) + a ab B (b) + e abc a bd B d E c - H ab E b j , (30) 

or 

B-Vx(txB) - AV 2 B = a ■ B + X ■ B - V x (a ■ E) + a ■ B + (a ■ V) x E - H ■ s] , (31) 

respectively, where H ab = (cmla) ab . The terms on the rhs describe the effect of the curvature, and given a GW they will 
drive the magnetic field evolution. We see that the finite electric conductivity gives rise to new couplings between GWs and 
magnetic fields. 
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4.1 Gravitational waves in a constant magnetic field 

The case of a GW passing through an initially constant magnetic field will help illustrate these equations. We treat all 
disturbances VB in the magnetic field as well as the gravitational wave variables as first order perturbations, and we therefore 
neglect all terms in the equations which are a product of derivatives of the magnetic field with GW terms (a ab ), and terms 
0([VB] 2 ). We may also neglect disturbances in the energy density of the fluid. With these approximations, we find that the 
momentum conservation equation l|16|l becomes 

Sv a = e abc j b B c , (32) 

while the current is given by Eq. 0: 

j a — (j,q 1 curlB a . (33) 

The dynamo equation 1301 simplifies to 

B a - e abc e cde B e D b v d - AD 2 B a = 2\& ab B b + o ab B b . (34) 
We may decouple Eqs. 132H - 1341 by taking the time derivative of Eq. 1341 : 

B a - AD 2 B a + -^B b [—B a D 2 B b + B c (D a D b B c - D b D c B a )] = 1\b ab B b + b ab B b (35) 

where (neglecting backreaction from the electromagnetic field) the shear perturbation is determined by the source-free wave 
equation 

■6 ab - D V ft = 0. (36) 

We split the perturbed magnetic field parallel and perpendicular to the 'background' magnetic field according to 

B a = B [(1 + 3S)e a + 38 a \ , 38 a = N ab B b /B , and 38 = (e a B a - 1) /B . (37) 

where N ab — h ab — e a e b and e a is a spacelike unit vector: e a e a = 1 fsee lciarkson fc Barrett] l|2003l) for a detailed discussion of 
this split). Here, Bo is a constant, denoting the magnitude of the magnetic field when the GW is zero (in which case B a — Boe a ). 
The perturbations in the static magnetic field are given by 38, parallel to the background field, and 33 a perpendicular to it. 
We neglect all products of 38 and 3§ a , and terms of the form 38a (note that Ve = 0(a)), since these are of higher order. 
The shear perturbation can be decomposed according to 

oab = E a6 + 2E (a e b) + (e a e b - \N ab ) E (38) 

by introducing the gravitational wave variables 

E = a ab e a e b , E a = N ab a bc e c , and E ab = (N( a °N b f ~ ^N ab N cd ^ a cd , (39) 

and the operators 

Dm = e a D a = — and = N ab B b , (40) 
dz 

allow us to split the spatial derivative along, and perpendicular to, the unperturbed magnetic field. 1 

Using these definitions, and the frame choice 2 e a — = D™ e a , we may split Eq. 135H into a longitudinal scalar equation 
along e a : 

J - AD 2 J - C 2 aD 2 3S = 2AE + E, (41) 

where Ca = (eoBo / '^) 1 ' 2 is the Alfven velocity of the MHD plasma, and a transverse vector equation perpendicular to e a : 

m a - XD 2 3S a - C\ {p\m a - DlD||^) = 2AE™ + E a . (42) 

We also have the D a B a = constraint, which takes the simple form T>^3§ + D^38 a = using our new variables. Neglecting 
backreaction, gravitational waves propagating along the background magnetic field will not affect the magnetic field evolution, 
due to the condition of vanishing divergence T> a <J ab = = e a a ab . This can also be seen from Eq. 138H . where (if e a represents 
the propagation direction of the gravitational wave) Y> ab is the only non- vanishing part of the decomposition of a ab , thus 
removing the coupling to the magnetic field via E and E a . On the other hand, including the backreaction from the magnetic 
field on the gravitational wave will in general make the gravitational wave non-transverse, by the introduction of small scale 
tidal forces as well as large scale curvature effects, but this issue is left for fu ture research. 

Next we introduce spatial harmonics for the perturbed variables (see Iciarkson fc Barrettl J2003)). Define the scalar 
harmonics 

D 2 ±Q( k± ) = -k 2 ±Q( k± ), (43) 

Note that t he definition of D° for ten sors has projection with N ab on each free index IClarkson fc Barret dEooah . 
2 We refer to lClarkson fc Barret J 120031) for a full discussion of frame choice vs. true GW degrees of freedom. 
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where k± is the harmonic number for the part of the fields perpendicular to e a . We can then define vector harmonics of even 
and odd parity, 

Q(k ± ) = (k±r 1 D° L Q (k±) Q\ k±) = (k x y 1 e a b D b ± Q ikx) , (44) 

respectively, where e a b = £abce c = —tba is the volume element on two-surfaces orthogonal to e a . These vector harmonics obey 
DlQ a (k±) = -k 2 ± Q a {k±) (similarly for Q a (k±) ). 

Using these definitions, we may expand the magnetic field variables as 

& = £^x)«(kL). (45) 
k ± 

® a = E^x)Q?^)+^L)Q(*x). (46) 

and similarly for the GW shear variables 

£ = ^2^(k ± )Q(k x ), and E a = ^ ^%k ± )Q(k ± ) + ^{k ± )Q(h ± )- (47) 

Using these harmonic expansions, Eq. 1 1H can now be written as a partial differential equation for each kj_, using ' = d z : 

23 s - C\®1 - + k\\3Bs + k\C\SS s = 2AE S + E Sl (48) 

while the even part of 1 121 may be derived from k±3?v = 2S'g (the div-B constraint), and is in fact the same equation, but 
with S replaced by V; these give the dynamics of the 'even sector', while for the odd sector we have 

JV - C\^v - AJv + A''" A. A + k\c\&9 v = 2AE V + t v . (49) 

The GW variables obey E - DjjE s + kj_Es = 0, and similarly for E v and Ey; D b <r ab = implies k ± T, v = DyEs- We 
have dropped k± subscripts, with the understanding that these equations apply for each value of k±. We may now introduce 
harmonics for <9 Z and dt in the usual way if desired. 



4.2 Integration 

As a means of understanding the effect of conductivity on a magnetic field in a GW spacetime we shall consider the case of 
a Gaussian pulse of GW travelling in the positive z-direction through a static magnetic field such that at t — 0, 

E^ fc) = Ee- z2/L2 cos(^), (50) 

E^ fe) = Ee~ z2/L2 Lsin(w«) + 2|-cos(w0)] , (51) 

where we set E = — 1 without loss of generality (alternatively, this could be Ey because the decoupled equation for 3$v is 
identical to that of 2§v)- This means that SS is measured in units of the GW amplitude in the plots below. We will now 
integrate Eq: 1481 for a variety of situations assuming that at t — 0, 3§s = 2§v = - i.e., the interaction switches on at t = 0. 
We take u — it/3, k± = (so 3§s ~ 0) and L = 8, below, all in units where c = 1. 
The qualitative features of finite conductivity are shown in figures and [5] 

In Fig. we consider the case of Ca = 1, that is waves in a very tenuous plasma, giving rise to resonant magnetic field 
amplification. The figures shown apply equally to SSy or 8$y. In the case of infinite conductivity, shown on the left, we have 
a linearly growing magnetic field pulse, which is continually sourced by the GW as it travels through the background field. 
When conductivity is present, however, this amplification is inhibited in that the induced magnetic field can only grow to a 
certain strength (compare the scales of the two plots) before a saturation amplitude is met, at an amplitude considerably less 
than if the magnetic diffusivity is zero. 

If the magnetic field is strong, or the fluid dense, then the Alfven velocity can be small, resulting in weak non-resonant 
amplification of GW. The plot on the left of Fig. shows the modes with Ca = 0.5 when A = 0. Dispersion shows up here 
as a trail of small amplitude waves, trailing behind the initial pulse, due to the mismatch in the GW and plasma dispersion 
relations. The plot on the right demonstrates what happens when the conductivity is finite. As in the resonant case the 
forward driven pulse is damped and can only reach a certain amplitude, which is considerably less than when Ca = 1- The 
trailing part of the wave has nothing to sustain it, so diffuses and rapidly decays away. 



5 DISCUSSION 

Finite conductivity can alter the behaviour of conducting fluids, plasmas and electromagnetic fields in a vast variety of ways. 
When spacetime curvature may be neglected these effects are well understood in a significant number of physical situations, 
ranging from the early universe, to astrophysical and laboratory plasmas. When spacetime curvature effects are present and 
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Figure 1. The effect of A on long wavelength modes, with Ca = 1- The plot on the left has A = showing linear growth of the magnetic 
field as the GW passes through the background field. On the right, we have the same situation but with finite conductivity, A = 0.05. 
While there is some early amplification of the magnetic field, this quickly saturates, with late time behaviour just a brief disturbance as 
the GW passes. Note the difference of the vertical scale between the two plots. 




Figure 2. The effect of A on slow Alfven waves, with Ca = 0.5. The plot on the left has A = showing little growth of the magnetic 
field as the GW passes through the background field, trailed by a slow wave. The dispersion takes the form of a slower wave trailing 
the oscillations induced by the GW. On the right, we have the same scenario but with finite conductivity, A = 0.05. The trailing Alfven 
waves are dissipated very quickly. The main pulse, as in the previous plot, is from a direct forcing from the GW. 



significant, however, very little is known about the role played by finite conductivity in astrophysics and cosmology. We present 
here, for the first time, the dynamo equation in its full generality in an arbitrary curved spacetime. 

The MHD approximation, describing low frequency charged fluid phenomena, allows us to neglect the displacement 
current in Ampere's Law, from which we may derive a diffusion equation for the magnetic field - the dynamo equation Hill 
or il l 2> . For easy co mparison with non-ge neral relativistic results we have utilised the intuitive 1+3 covariant approach to 
relativistic analysis ijEllis fc van Elstlll998l) . This approach splits spacetime into space and time in a covariant way by use of 
a timelike vector field, which then also plays the role of the convective derivative, but now on an arbitrary curved manifold. 
Spacetime dynamics is then covariantly described by the use of invariantly defined scalars, 3-vectors, and PSTF tensors, made 
up from the Riemann curvature tensor, and the dynamical quantities associated with u a . These quantities feed into Maxwell's 
equations in a non-trivial way, making the resultant electromagnetic field experience 'gravitational currents' of a very different 
nature from flat space physics. 

In the relativistic dynamo equation we have derived, these gravitational currents come in a huge variety of terms which 
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are very difficult to analyse in general. We have presented in some simplifying assumptions which are often used in spacetime 
modelling, in order to better gain an intuition about the kinds of effects conductivity coupled to gravity can produce. 

As a specific example of analysing the relativistic dynamo equation, we considered in detail the case of a plane GW 
interacting with a static magnetic field. We believe this may play an significant role around compact objects, where magnetic 
fields and gravity waves can interact strongly. By introducing covariant vector harmonics adapted to the magnetic field we 
reduced the dynamo equation to two fairly simple coupled PDEs with a constraint, and one decoupled PDE, with the space 
dimension along the background field. These PDEs for the induced magnetic field are sourced by the GW in an intuitive 
fashion, and may be analysed by standard means. We integrated these equations numerically for a Gaussian pulse of radiation 
passing through the background static magnetic field; the summary of this integration is presented in the figures. In the case 
of infinite conductivity, there is linear amplification of the magnetic field which mimics the GW; this is followed by a trail of 
oscillations resulting from waves in the magnetic field travelling with the Alfven speed. When magnetic diffusivity is present, 
these trailing modes diffuse quickly, and the overall amplification of the magnetic field is inhibited; a ripple in the magnetic 
field is forced as the GW passes through. 
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